Interaction-induced decoherence in non-Hermitian quantum walks of ultracold Bosons 
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We study the decoherence caused by particle interaction for a conceptually simple model, a 
quantum walk on a bipartite one-dimensional lattice with decay from every second site. The corre- 
sponding non-interacting (linear) system has been shown to have a topological transition described 
by the average displacement before decay. Here we use this topological quantity to distinguish 
coherent quantum dynamics from incoherent classical dynamics caused by a breaking of the trans- 
lational symmetry. We furthermore analyze the behavior by means of a rate equation providing a 
quantitative description of the incoherent nonlinear dynamics. 

PACS numbers: 03.65Vf, 03.65 Yz, 03.75.Gg, 05.60.Gg, 64.70Tg 



The quantum-to-classical transition often denoted as 
decoherence is of increasing recent interest and consid- 
erable progress has been achieved 0,0. Open many- 
particle systems, in particular ultra-cold quantum gases 
and Bose- Einstein condensates, are well-suited for inves- 
tigating decoherence processes in nonequilibrium quan- 
tum dynamics. Recent examples of theoretical and 
experimental studies of such systems include quantum 
walks in optical lattices and the decay dynamics in 
open, non-hermitian systems [ 6rjl3l . some of which ex- 
plicitly address decoherence [llHl3j . 

In this brief report, we analyze the mechanism of deco- 
herence in the nonequilibrium dynamics of an interacting 
Bose gas by means of a topological quantity introduced in 
[IH ■ Topological properties have also been of interest in 
other recent investigations with cold atoms concerning, 
e.g., quantum Hall effects [IE El, Berry phases fl7l Hq 
or dissipative quantum wires [19(. We will demonstrate 
an interaction induced decoherence, i.e. a quantum-to- 
classical transition for a conceptually simple model sys- 
tem. In particular, we study the quantum walk of ultra- 
cold bosons in a deep bipartite ID optical lattice with two 
different tunneling rates and decay from every second site 
described in a non-hermitian ti ght- binding model con- 
sidered by Rudner and Levitov jl4| in their analysis of 
a topological transition. Here we generalize and study 
the non-hermitian discrete Gross-Pitaevskii (or nonlin- 
ear Schrddinger) equations 
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where the nonlinear terms g\a m \ 2 , g\b m \ 2 model the inter- 
action between the particles. (Note that such an equation 
introduced hcuristically in studies of ultra-cold bosonic 
gases or optical devices poj can be derived rigorously as 
a mean-field approximation from a multi-particle Bosc- 



Hubbard model with dissipation in the limit of high par- 
ticle numbers in the system fl2j . In a different approach, 
however, a mean- field limit with a renormalized interac- 
tion term appears [2l],|2^].) The system ([1]) is illustrated 
graphically in Fig.Q] Initially all particles occupy a single 
non-decaying site which corresponds to the initial condi- 
tions a m (0) = S m o, 6 m (0) = at time t = 0. 




FIG. 1: (Color online) Bipartite ID lattice with couplings v 
and v' and decay from every second site with decay rate 7. 

In the noninteracting (linear) case, ([!]) with g = 0, 
Rudner and Levitov [14| analyzed the topological quan- 
tity 



pOO 

Am = y^TOP m , P m = / 7|6„ 

Jo 



{t)\ 2 dt (2) 



describing the average displacement of the particles be- 
fore their decay. By means of a Fourier transformation to 
momentum space they showed that Am is equal to the 
winding number of the relative phase between compo- 
nents of the Bloch wave function which leads to a quan- 
tization of Am as a function of the ratio v/(v + v') of 
the tunneling coefficients: Am = 1 for vj (y + v ') < 1/2 
whereas Am = for v/(v + v') > 1/2. In [14| the quan- 
tization of Am was shown to be robust against intersub- 
lattice dephasing and classical noise on the energy levels 
e a , eb by means of numerical simulations. This quantized 
behavior was compared to the continuous dependence of 
Am on v and v' in the case of incoherent tunneling. 

The latter can be obtained by assuming a randomly 
fluctuating phase between the sites so that the rates for 
incoherent hopping are given by perturbation theory to 



the lowest order as T oc v 2 , V 



[gj . Initially only 



the site clq is occupied. After the very first hopping pro- 
cess to one of the decaying neighbour sites bo or b\ corre- 
sponding to m = and m = 1 respectively, the incoher- 
* Electronic address: rapedius@physik.uni-kl.de, korsch@physik.uni-kl.dent dynamics between the decaying sites b m is symmetric 



2 



with regard to transport to the left and right so that the 
contributions of all subsequent hoppings to the displace- 
ment (|2|) cancel out. Since the m = term does not 
contribute, the displacement Am is given by the proba- 
bility of an initial hopping to site b\ with m = 1 which 
is equal to the corresponding relative hopping coefficient, 
i.e. 

Am incohcrcnt = r'/rr + r') = v l2 /(v 2 + v' 2 ) . (3) 

Figure [2] compares the the quantized behavior of Am 
for coherent tunneling (dashed line) with the continuous 
curve §3§ (solid line) for incoherent tunneling. 

In order to investigate the influence of a finite mean- 
field interaction g ^ we integrate the system (TT]) nu- 
merically on a finite lattice with periodic boundaries and 
parameters e a = 65 = 0. In all numerical calculations 
scaled units with K = 1 are used throughout this paper. 
The results are also displayed in Fig. [2] Even in the lin- 
ear case .9 = the numerical results for Am deviate from 
the perfectly quantized behavior of the idealized infinitely 
extended system around v/(v + v') » 0.5 due to effects 
of finite size and finite integration time. These effects 
were discussed in detail in [lj] and are not significantly 
altered by the mean-field interaction. For increasing val- 
ues of the interaction parameter g, the values of Am 
deviate more and more from the quantized behavior of 
the linear system, approaching the incoherent tunneling 
curve for higher values of g. Thus the mean-field interac- 
tion, unlike other kinds of disturbances mentioned above, 
induces decoherence in the system. 

The reason for this decoherence lies in the fact that 
the nonlinear interaction breaks the translational sym- 
metry within the non-decaying sublattice as will become 
clear in the following. To support our argument, we will 
demonstrate that, instead of considering interactions, the 
coherence observed in the linear system can also be dis- 
turbed by means of a much simpler symmetry breaking 
mechanism, namely by adding a constant energy shift to 
the initial site ao . Figure [3] shows the numerically cal- 
culated displacement Am as a function of v/(v + v') for 
the system {T]) with g = and an additional energy shift 
77 of the initial site energy, i.e. we make the replacement 
( a — > £ a + vfimo in © • We observe that the values of Am 
more and more approach the incoherent hopping curve 
for increasing values of the energy shift 77. By means of a 
constant energy shift of a single site we are thus able to 
mimic the decoherence inducing effect of the mean-field 
interaction, thereby identifying the breaking of the trans- 
lational symmetry within the non-decaying sublattice as 
the main cause of decoherence. 

To gain further insight we analyze the time-dependence 
of the correlations a*(t)a m (t), £>*(t)6 m (t) and 6*(t)a m (t) 
between the sites within and between the two sublattices 
which is most conveniently done by reexpressing our sys- 
tem by means of a density matrix. At first we concentrate 
on the simpler case of a non-interacting system with an 
additional shift of the initial site. For convenience, we 
first rewrite the system |T]) with g = and an additional 



0.5 




v/(v+v'; 



FIG. 2: (Color online) Displacement Am as a function of 
the ratio v/(y + v') of the tunneling coefficients. '- -': linear 
(3 — 0) infinitely extended system, ' — ' incoherent hopping 
according to @. The symbols represent numerical results for 
a finite system with 46 lattice sites and periodic boundary 
conditions: 'X': 3 = 0, 'A': 3 = 0.2, V: g = 0.5, g = 1 
,'0': 3 = 4. The parameters in the numerical simulation are 
7 = 2, v' = 0.5 with < v < 0.5 for v/(v + v') < 0.5 and 
v = 0.5 with < v < 0.5 for v/(v + v') > 0.5, e a = = e b . 
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FIG. 3: (Color online) As Fig. [2] but for the linear system 
(3 = 0) with an energy offset 77 of the initial site. 'X': rj = 0, 
'A': 77 = 0.05, V: 77 = 0.1, '0': V = 0.6. 



energy shift 77 of the site ao in the compact form 

\hc a = -J' a c a+ i - J a c a -\ + E a c a (4) 

where 



a m=a /2i a even E _ ( e a + 7jS a0 , a even 
&m=(a+i)/2i a odd ' a \e h -ij/2, a odd 

(5) 



and 



T _ J v/2, a even , _ J v'/2, a even , . 

Ja ~ \ v'/2, a odd ' Ja ~ \ v/2, a odd ' (b) 

For the dynamics of the density matrix elements p a p = 



3 



c* C/3 Eq. ((U) then yields 

\hp a p = (E/3 — E^Pap + J' a Pa + l,l3 + JaPa-1,/3 

-J'pPa,f} + l ~ JpPa,/3-l ■ (7) 

The initial conditions at t = read p a p(0) = S a p5 a Q. 
As above we assume e a = = e& in the following. We 
consider the dynamics of the central diagonal element of 
the density matrix 

ihpoo = JoOio - Poi) + Jo{p-i,o - Po,-i) (8) 

which depends on the nearest off-diagonal elements. As 
an example we have a closer look at the element 

ifi-Poi = (Ei - Eq)pqi + JqPh + JoP-i.i - J1P02 - JiPoo ■ 

(9) 

If we neglect the off-diagonal elements between two odd 
(i.e. decaying) sites we arrive at 



ihp i « (-Ei - Eq)p i + Jq(pii 
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J0P02 (10) 



where we have used J\ = J' . The dynamics of the 
off-diagonal clement P02 must be examined more closely. 
Since our initial state is localized at site wc can neglect 
decaying off-diagonal elements not involving site 0. If we 
additionally neglect correlations over a distance of three 
or more sites Eq. ([7]) yields ihp 02 ~ {E 2 - Eq)p 02 - Jo Poi- 
Assuming a slow time-dependence of poi compared to 



h/(E 2 



we obtain P02 ~ Poi(l — exp(— i(E 2 



E^)t/h))J /(E 2 -E^. For \E 2 -E^\ = |t?| » J the term 
J0P02 in (HOI) can thus be neglected. Physically this cor- 
responds to a suppression of resonant tunneling between 
the sites and 1 due to the large difference \E 2 — E^ \ = \rj\ 
of the on-site energies. Following the reasoning in [24| we 
assume a slow time-dependence of poo and pn compared 
to the real part of h/(Ei — Eq), which is justified for 
\Re(Ex - E%)\ = \rj\ > J , to obtain 



Poi 



J(i Pn ~ Poo 
h 77 + 17/2 
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(11) 



with Ei — E^ = —i] — vy/2 and, analogously po,-i = 
{J /K){p- h -i - p o)/(r? + i7/2) - P* i,o- Thus the dy- 
namics of the central diagonal matrix element ([5]) be- 



comes 



Poo 



-(r + r')p 00 + rp_i._! + r'pn 
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where 



r = 



(JoA) 2 7 _ (v/h) 2 j 



rf + 7 2 /4 Arf + 7 2 



r' = 



(J„A) 2 7 (v'/hfj 



V 2 + 7 2 /4 4?y 2 + 7 2 ' 
(13) 

Note the quadratic dependencies T oc v 2 , V oc v' 2 already 
stated before Eq. ©. For a sufficiently strong decay 
coefficient 7 (compared to the tunneling coefficients v and 
v') one can assume that the correlations between sites 
further away from site 0, which vanish at t = 0, do not 
build up significantly in the course of the decay process. 



As an approximation we may thus assume an incoherent 
dynamics as described by (fT2"j) for all sites with some local 
site-dependent transition rates r a and T' a yielding a rate 
equation 



p aa ~ -(r + r')po, Q 

^ f r Q p Q _i, Q _i + T' a p a+lt0t+ i, a even 

1 -JPaa + ^' a Pa-l,a-l + ^aPa+l,a+l, a odd. 



(14) 



However, due to the decay in the system the quantities 
that we are interested in, namely Am and the occupation 
of the central site poo (cf. below) are not sensitive to the 
exact dynamics in the outer sites whose main effect in this 
context is the suppression of boundary effects. Thus for 
simplicity we make the approximation of constant tran- 
sition rates r Q — T and r' Q = T' in the following. Within 
this approximation, equation (|14p can be solved in closed 
form in Fourier space with the result ([5]) for Am (see Ap- 
pendix). 
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FIG. 4: (Color online) (Color online) Dynamics of the linear 
system (g = 0) with an energy offset 77 = 0.6 of the initial 
site a = for the parameters v = 0.25, v' — 0.5, 7 = 2 and 
e a = = £6. The numerical integration of the full system 
(J4j) (poo:solid line, Am t : bold solid line) is compared with the 
dynamics according to the rate equations Q14[l (poo^dashed 
dotted line, Amt: dashed line) . 



In order to compare the dynamics obtained from the 
rate equation (|T4"|) with the dynamics of the full system 
we define the time-dependent displacement Am t , which 
is obtained if the time integration in ([3]) is only performed 
up to a finite time t such that Am = lim^oo Am t . An 
example of the dynamics of Amt for ?; = 0.6 is shown 
in Fig. 2] together with the corresponding decay of the 



central site occupation poo 



l a o| 



While poo decays 



exponentially, the displacement Am t increases until it 
reaches its final value, given approximately by for 
long times. For both quantities the numerically exact 
calculation is reasonably well approximated by the rate 
equation result. 

Let us now return to the original, nonlinear problem 
with a finite interaction g ^ 0. The role of the energy 
offset 77 is now played by the mean-field interaction. For a 
sufficiently strong decay coefficient 7, the instantaneous 
energy offset between the central site and its neighbors 
is approximately given by the local mean-field interac- 
tion term at site 0. To obtain an approximate descrip- 
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a better understanding of the transition from coherent 
quantum dynamics to incoherent classical dynamics in 
mesoscopic systems. 

We thank Eva-Maria Graefe and Dirk Witthaut for 
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FIG. 5: (Color online) As Fig.Ubut for the nonlinear system 
with g — 4 (upper panel) and g — 0.5 (lower panel). 



Changing the notation to pm = Paa for a even (+) or 
odd (-) Eq. (TT4")) with constant decay rates r a = V and 
r' = F reads 



(+) 



-rt + r'77 (+) 



m+l 
( + ) 



(15) 



with r = T + T' and T' Q = T + 7. As in [14J we switch 
to the momentum representation 



,(±) 
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X>- ifen p£>, (16) 



tion of the nonlinear system dynamics we thus use the 
rate equation (fhfj) with the replacement 77 = gpao(t)- 
Even though the hopping rates T and V are now time- 
dependent, their ratio remains constant on this level of 
approximation such that the final displacement Am is 
still well approximated by The upper panel of Fig. [5] 
demonstrates that for g = 4, corresponding to the inco- 
herent regime (cf. Fig. [5J , the modifications of the system 
dynamics due to the nonlincarity are well described by 
the effective rate equation, both for the time-dependent 
displacement Am* and the central site occupation pqq. 
The latter now shows a nonexponential decay behavior 
typical of open nonlinear systems ■ ^ u ^ ne l° n g time 
limit, the exponential decay of the linear system is re- 
covered as the influence of the nonlinear interaction term 
<7Poo(t) becomes negligible. For comparison the lower 
panel of Fig. [5] shows the dynamics for a smaller interac- 
tion g = 0.5 corresponding to an intermediate, still par- 
tially coherent regime (cf. Fig. [2} . While the decay of the 
central site occupation is still described reasonably well 
by the rate equation approach we find clear deviations 
for the time-dependent displacement Arrit as expected. 

In summary, the interaction-induced decoherence of 
cold bosons spreading in a bipartite optical lattice with 
decay was analyzed by means of a topological quantity. 
Using numerical and analytical methods, the breaking of 
translational symmetry within the non-decaying sublat- 
tice and the resulting suppression of resonant tunneling 
was identified as the cause of decoherence. For the regime 
of strong decay and strong interaction, a rate equation 
providing a quantitative description of the system's in- 
coherent nonlinear dynamics was derived. The authors 
hope that the present model study can contribute to 



where the integration extends over the Brillouin zone 
— 7r < k < 7r. The resulting equations for the momen- 
tum distributions 



q k - -log* +-lfe9fc 



*(-) - 



(+) 



(17) 



with Tfc = r + r' e can be solved immediately and a 
solution for the initial conditions q^\o) = 1 , g^ ^O) = 
(the translation of Pm\o) = S m o , Pm\o) = 0) is 



V = 



+ ^ r ((A + +r[ ) )e A + t -(A_+r^)e A - t ) 

(18) 



e A + t_ e A_t 



with 



a± = -(r + r^)/2± V7 2 /4 + |r fc | 2 . (19) 

The momentum representation of (J2J) is given by 

p 00 poo 

Am = 7/ dtY,mp(-){t)= 1 / dtid k qir\t)\ k=0 

poo 

= Hd k Qk\ k „ with Q k = dtq[r\t). (20) 

Jo 

Integration of the solution q k (t) in (|T8| yields 



Qk = 



1 k 



and Am 



r 7 



r + r' 



(21) 
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